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Data from 1 1 years of continuous spectroscopic observations of the active RS CVn-type binary star EI Eridani - gained at 
NSO/McMath-Pierce, KPNO/Coude Feed and during the MUSICOS 98 campaign - were used to obtain 34 Doppler maps 
in three spectroscopic lines for 32 epochs, 28 of which are independent of each other. Various parameters are extracted from 
our Doppler maps: average temperature, fractional spottedness, and longitudinal and latitudinal spot-occurrence functions. 
We find that none of these parameters show a distinct variation nor a correlation with the proposed activity cycle as seen 
from photometric long-term observations. This suggests that the photometric brightness cycle may not necessarily be due 
to just a cool spot cycle. The general morphology of the spot pattern remains persistent over the whole period of 1 1 years. 
A large cap-like polar spot was recovered from all our images. A high degree of variable activity was noticed near latitudes 
of ^60-70° where the appendages of the polar spot emerged and dissolved. 
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1 Introduction 

Our Sun exhibits cyclic behaviour which can be seen in its 
spectral output as well as in its spot coverage ( [Willson & 
[Mordvinov 20Q3| ). Apart from the commonly known 11- 
years Sunspot cycle, our Sun shows an 80 years long so- 
called Gleissberg cycle and an 200-300 years long pseudo- 
cycle (Wolf, Sporer, Maunder). Even longer cycles or vari- 
ations are seen, and there is evidence that periods of cyclic 
behaviour alternate with times of constant, very low activity 
(e.g. Maunder minimum). Solar-type G and K stars can ex- 
hibit the same solar-like cycles in their photospheric output 
as originally discovered from their chromospheric Ca II H 
& K emission (| Baliunas & SQon|p^95l [Gray et al.|p^ ). 
However, a photosphere-chromosphere correlation cannot 



be easily reproduced for evolved G and K stars ( Choi et al. 
1995[ |Frasca et"aL]|2005| ). Long-term changes in the mean 



brightness of RS CVn stars are therefore not a straightfor- 
ward tool for investigating spot cycles similar to the Sun's 
1 1-year cycle. From what we know from our Sun, we would 
expect an activity cycle from Ca II H & K to be accompa- 
nied by a spot cycle. 

Several groups began long-term Doppler-imaging stud- 
ies on a few selected active stars (that usually consist of 
one or a few images per year, i.e. per observing season): 
Korhonen et"aL| ( [20071 Com, 24 maps from 1993 to 
2003), IVogt eFaL] ([T9991 HR 1099 = V711 Tau, 23 maps 
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from 1981 to 1992),'Berdyugi na et al.| ( [T998| II Peg, 9 maps 
from 1992 to 1996; and 2001 , LQ Hya, 9 maps from 1993 
to 1999), Olah et al. ( 2002b, UZ Lib, 8 maps from 1994 to 
2000), Marsden et al. (j2007| IM Peg, 31 maps from 2003 
to 2005), Donati et al!l pOQ3| AB Dor, LQ Hya, HR 1099) 
and others. All of them are active stars and most of them are 
close binaries with their rotation period usually spun up by 
tidal forces. However, so far only one star, HR 1099 (| Vogt 
|& Hatzes||1996||Strassmeier & Bart"usl|200Q| , has been ob- 
served over a sufficiently long time span to cover a complete 
activity cycle with Doppler maps. 

EI Eridani = HD 26337 (G5 IV, P^ot = 1-947 days, V 
= 7.1 mag, vsini = 51 kms~^, SBl) was one of the first 
stars being Doppler imaged. Its spot distribution has been 
monitored since 1984. Long-term photometric observations 
revealed large brightness variations which seem to be cyclic 
and were first supposed to resemble a solar-like 11 years 
cycle but appear to be more complex as observations con- 
tinue. Within the first 16 years of photospheric observations, 
EI Eri gave the impression of a cycle of similar length to the 
solar cycle: [Strassmeier et al.| ( |1997] ) calculated a length of 
11 ± 1 years. Observations in the following years did not 
confirm this cycle as they did not meet the anticipated de- 
cline of brightness. [Berdy ugina & Tuominen| (p^998 ) esti- 
mated a 9-year periodicity from the positions of active re- 
gions. They did not discuss a possible change in the level 
of spottedness itself, though. [Olah et aL| ( |200 0|) favoured a 
length of 16.2 years on a time base of 18.5 years and quote a 
significant remaining period of 2.4 years. Adding four more 
years of data, the preferred cycle length was back to 12.2 



years (see Olah & Strassmeier 2002). Most recent calcu 
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lations, comprising an observed time interval of 28 years, 
state a long cycle of about 14 years with variable ampli- 
tude and two short cycles of about 2.9-3. 1 years and 4. 1-4.9 
years ( Olah et al.|2009| ). 

The present program intends to review all available 
Doppler maps of EI Eri and to possibly link any long-term 
evolution of the spot distribution to the photometric varia- 
tions. 

2 Doppler imaging of EI Eridani 

With its large rotational velocity and an intermediate incli- 
nation, EI Eri is, as |Fekel et al.| ( p^87) ) already noted, an 
ideal candidate for Doppler imaging (hereafter "DI"). Con- 
sequently, EI Eri has been a prime target since the first ap- 
plication of this technique to spotted late-type stars in 1982. 
Doppler maps can be found in |Strassmei"er| ([199Q| epoch 



Spectrum of EI Eri 



1987), [Stessmeier et al.] ( p^9T| epoch 1988; investigation 
of different DI techniques), and for the 1984-87 period, in 
Hatzes & Vogt| ( |1992| ). |Q'Neal et al.| ( |1996| ) determined spot 
covering factors between 16% and 37%. 

The data used for producing the maps presented in this 
paper were obtained at NSO McMath-Pierce during seven 
years of our long term synoptic observations (1988-1995; 
19 maps) and one dedicated visitor observing run in 1996 
(three independent maps), several dedicated KPNO/Coude 
feed visitor observing runs in 1995, 1996 and 1997 (the lat- 
ter with three independent maps) and participation in the 
MUSICOS 1998 observing campaign (three independent 
maps). The total covered time period streches from 1988 
to 1998 (eleven observing seasons). A detailed description 
of the NSO and KPNO observations and the obtained data 



was given in [Washuettl et al.| ( |200 8 , hereafter paper I). The 
observations and data from the MUSICOS 98 campaign will 
be presented in a forthcoming paper ( Washuettl et a l.|2009| 
which will also address the question of differential rotation. 
Doppler imaging suggests abundances of —6.3 dex (Ca) and 
—5.5 dex which are smaller than solar by —0.6 dex (Ca) and 
— 1.1 dex (see below). 

Knowledge of the precise stellar parameters is essen- 
tial for obtaining good-quality Doppler maps. Parameters 
like temperature, rotational velocity, luminosity, radius, or- 
bital inclination, mass, period and orbital parameters were 
already investigated in paper I. Additionally, several specific 
Doppler-imaging parameters were improved by producing 
a series of test reconstructions: stellar inclination, differen- 
tial rotation, rotational velocity, gravity, temperature, abun- 
dances, macro and micro turbulence, line blends and transi- 
tion probabilities. Nomalized distributions for these pa- 
rameters which allow to find or improve the parameter value 



were presented in , Washuettl ( 2004|. For a list of Dl-relevant 
parameters, see Table [T] Further as trophy sical parameters 
are listed in paper I, Table 5. 

For DI, it is important to obtain high signal-to-noise ra- 
tios (S/N > 200) at moderate to high spectral resolution 
(A/AA ^ 30-40 000). The wavelength region around 6420 




Fig. 1 Typical spectrum of EI Eri with the four mapping 
lines A6393, A6411, A6430 and A6439. 



Table 1 Values used as input parameters for the Doppler 
imaging code TempMap. 



Parameter 


Value 


-Pphot — Prot 


1.9472324 


To, phot =To,rot 


2448054.7109 


7 


21.64 kms"^ 


Kl 


26.83 kms"^ 


e (eccentricity) 


(adopted) 


^phot 


5500 K 


T 

-L max 


6000 K 




51.0 kms"^ 


Inclination i 


56.0° 




3.5 


Micro turbulence ^ 


2.0 km s"^ 


Macro turbulence Cr = Ct 


4.0 km s"^ 


Regularisation 


Maximum Entropy 


Weight of phot, data 


0.1 


log [Ca] abundance 


—6.3 (0.6 dex below solar) 


log[Fe] abundance 


—5.5 (1.1 dex below solar) 



A provides up to four relatively unblended lines usable for 
DI: Cal 6439, Fel 6430, Eel 6411, and Eel 6393 and was 
therefore chosen for the single-order spectrographs avail- 
able at the McMath-Pierce and Coude feed telescope. An 
example spectrum is shown in Fig.[T] 

EIEri's unfortunate rotational period of 1.947 days is 
almost an integer multiple of the day/night cycle. As a con- 
sequence, one ideally needs 20 nights of continuous obser- 
vations from a single observing site in order to achieve per- 
fect phase coverage. In practice, 14 consecutive nights are 
sufficient to give a good-quality Doppler image. No correc- 
tion for the contribution of the companion star had to be 
adopted as EI Eri is a single-lined spectroscopic binary. The 
time resolution of consecutive spectra was 2700 - 3600 s, 
corresponding to a spot motion of 6 - 8° in longitude on the 
central meridian of the stellar surface or 5.2 - 6.9kms~^ 
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Fig. 2 Comparison of all Doppler images from the Ca I 6439 line, consisting of pole-on view (left part), with dark, solid 
circles drawn at equally spaced intervals of 30° down to a latitude of -30°; and pseudo-Mercator projection (right part), 
from latitude -56° to +90°. The phases of the observations are marked by arrows around and below the maps, respectively. 
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(0. 1 1 - 0. 15 A) in the spectral line profile. This is still within 
one resolution element across the broadened line profile: 
With a typical resolving power of 24 000 - 36 000 (A/ AA), 
we achieve 8-12 resolution elements across the stellar 
disk, which corresponds to a velocity resolution of 8.3 - 
12.5 km (0.18 - 0.27 A) and a spatial resolution along 
the equator at the stellar meridian of approximately 10-14°. 
The orbital smearing accounts to 2.7 - 3.6kms~^ (0.06 - 
0.08 A). We therefore extended our integration limit to 60 
minutes. 

For the maps presented in this investigation, we apply 
the Doppler-imaging code TempMap by John Rice - as 
described by Rice et aL]([1989| ) and reviewed by Piskunov 
& Rice| ( |1993| ), |Rice| ( |1996| ) and, most recently, by Rice 
( |2002| ). All maps were plotted using the same temperature 



scale from 3600 to 5500 K. The spot temperature was de- 
termined by |Strassmei"er| ( [1990| ) using standardised V and 



R photometry and the Barnes-Evans relation (Barnes et al. 
1978) ). He gives a temperature difference (star minus spot) 
of 1860 ± 400 K and, with Tstar = 5460 K, yields an effec 



live spot temperature of 3600 ± 400 K. |Q'Neal et aLl ( [T996| 
derived spot temperatures by using the 7055 and 8860 A 
bands of the titanium oxide molecule and determined, as- 
suming a photospheric value of Tphot = 5600 K, a value 
of Tspot = 3700 ± 200 K which was confirmed by O^Neal] 



et al.| ( jl998) . In practice, the spot temperature is automati- 
cally chosen by TempMap and reaches down to 3600 K. 

All 30 independent maps for the A6439 line are shown 
in Fig. [2] Like most maps, it resembles a steady pat- 
tern: an asymmetric polar spot is escorted by several small 
low-latitude features. In some cases, spots slightly warmer 
(^200 K) than the photosphere occur. These are likely ar- 
tifacts which are produced by the Doppler imaging code 
due to, firstly, the use of differential photometry instead of 
absolutely calibrated photometry and, secondly, due to ex- 
ternal uncertainties in the spectra (see |Rice & Strassmeier| 
|2000| ). Lower metallicity values lessen the appearance of 
artifical hot spots. Therefore, we assume the high diver- 
gence of the Fe abundance from the solar case of —1.1 dex 
to be overdone and favour a more modest value of around 
—0.5 dex. This value is in better agreement with the evolu- 
tionary tracks by |Pietrinfemi et al. (2004) see paper I) and is 
also supported by Nords trom et al .| ( |2004| ). However, recom- 
puting the Fe Doppler maps with the more modest value of 
—0.5 dex does not alter the surface structure as such but only 
increases the temperature of the artificial hot spots. There- 
fore, we used the lower metallicity value of —1.1 dex for 
computation of the Fe maps in order to improve their com- 
parability. 



3 Analysis and parameterization of Doppler 
maps 

The spectroscopic data were inverted into a series of 
Doppler images spanning 11 years and amounted to a to- 
tal of about 100 images for up to 28 independent epochs (34 



images in A6439, 30 in A6430 and 25 in A641 1 and a few in 
A6393). Mean separation in time between two independent 
maps is 141 days, minimum and maximum separation is 3 
(overlapping) and 382 days, respectively. 
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Fig. 3 Equivalent widths for the three mapping lines at all 
epochs. The large open symbols are the mean values for the 
data used in the respective Doppler maps, the small symbols 
denote the individual spectra. The black lines denote linear 
fits for each mapping line. The corresponding values of the 
fits are listed in Tab. [2l 



Apart from tagging preferred spot locations, we extract 
several parameters from the Doppler maps and carry out sta- 
tistical evaluations to see if there are any trends in the tem- 
perature distribution or even correlations with the proposed 
activity cycle (see Fig. 2c in Qlah & Strassmeier|2002 and 
Fig.[7]in this paper). 



3.1 Line equivalent widths 

Fig. |3] shows the equivalent widths (EW) for each map- 
ping line (small symbols) as well as the average EW (big 
symbols) for each set of observations. A possible trend 
is the linear increase of the EW during the course of the 
eleven years of observation. This trend is seen with almost 
the same amount of 11.2 and 11.3% (2.4 and I. la) in the 
high-quality line A6439 and the mid-quality line A641 1, but 
only 4% (0.8a) in the poor-quality A6430 line (see Tab. [2]). 
If real, this could be suggestive of an increasing chromo- 
spheric contribution from a magnetic cycle with a period 
much longer than our observational coverage. 



3.2 Surface temperature 

Fig. |4] shows the progression of the average temperature for 
all Doppler maps from 1989 to 1999, distinguished for each 
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Average temperature of all Doppler maps 

1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 



Examples of temperature thresholds in Doppler maps 

A6439: 91/11 94/03 98/Ola 




Fig. 4 Temperature range for all Doppler maps. Filled cir- 
cles give the overall average of the respective map, bars de- 
note the minimum and maximum temperature. The solid 
lines are the low-latitude averages (-55° - 0°), the dashed 
lines the mid-latitude averages (0° - 45°), and the dotted 
lines the high-latitude averages (45° - 90°). 

Table 2 Increase of equivalent width, from Fig.js] 



Wavelength 


A6439 


A6430 


A6411 


Initial value (fit) [A] 


0.201 


0.172 


0.143 


Final value (fit) [A] 


0.223 


0.179 


0.159 


Increase value [A] 


0.022 


0.007 


0.016 


Increase percentage 


11.2% 


4.2% 


11.3% 


Standard dev. 


0.009 


0.009 


0.010 


Incr. value/Std. dev. 


2.4cr 


0.8cr 


1.7cr 



Table 3 Increase of fractional spottedness, from Fig. ^ 
Units are percentage points. 



Wavelength 


A6439 


A6430 


A6411 


Initial value (fit) 


11.2 


15.0 


16.6 


Final value (fit) 


12.6 


21.5 


18.9 


Increase value] 


1.47 


6.43 


2.24 


Increase percentage 


13.2% 


42.7 % 


13.5% 


Standard dev. 


2.7 


3.4 


1.8 


Incr. value/Std. dev. 


0.54cr 


1.9cr 


1.2cr 



mapping lin^ The filled circles are the overall temperature 
averages of the corresponding Doppler maps, the bars in the 
axis of ordinates denote the minimum and maximum tem- 
perature, the respective upper solid line is the temperature 
average of the low-latitude band, the dotted line that of the 
mid-latitude band and the lower solid line that of the high- 



^ Doppler maps of the A6393 Hne are omitted in the following inspec- 
tions as there are only four maps of A6393 available (one in 96/01 and three 
in 98/01). 
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Fig. 5 Temperature threshold definitions for three exam- 
ples of temperature distributions in our Doppler maps (Nov. 
1993, Mar. 1994, Jan. 1998). The arrows represent several 
threshold definitions for calculating the fractional spotted- 
ness in Fig.[6j 



latitude band. The low, mid, high-latitude bands run from 
-55° to 0°, from 0° to 45°, and from 45° to 90°, respec- 
tively (smaller latitudes are not seen as they are tilted out 
of sight due to the stellar inclination). Temperature varia- 
tions are noticed but not resembled in all spectral lines at 
the same time and therefore are likely not real. However, 
there seems to be a spot maximum (corresponding to a min- 
imum in temperature) in the high-latitude band in 1992 and 
again in around 1995, seen in A6439 and A6430. Overall, 
no pronounced correlation with the photometric variation is 
evident though. Systematically wrong input values for the 
line-profile computations can lead to spurious variations in 
the overall temperature, e.g. caused by artificial hot spots. In 
order to account for possible misfits, the latitude-dependent 
temperature averages were normalized by the respective to- 
tal average temperature. More clearly, we now perceive the 
temperature decrease in high latitudes in 1992 (in the A6439 
and A6430 line) which corresponds to an increase in the 
low-latitude band. We also see a clear overall negative cor- 
relation between temperature variations in the high and the 
low latitude bands. However, as in the pre-normalized case, 
all variations seen are not resembled in all spectral lines and 
therefore are likely not real. Above all, we do not perceive 
an apparent correlation between temperature and any other 
of the parameters derived from the Doppler maps. In par- 
ticular, no correlation with the long-term photometric vari 
ations as seen in Fig. [7] and, more clearly, in |01ah & Strass 
Fig. 2c), is evident. 



meter 



(2002 



3.3 Fractional spottedness 

We define a temperature threshold value and count the num- 
ber of pixels on the star's surface with temperatures below 
this threshold and relate it to the total number of surface pix- 
els. Several options for defining a threshold are presented in 
Fig. [5] which displays three examples of temperature distri- 
butions from the A6439 line. First, the overall temperature 
average is established, which usually sits at the lower edge 
of the photospheric temperature domain (5000-5500K). Its 
location is, however, so close to the edge that, as a second 
threshold option, a reduction of the temperature average by 
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Fig. 7 Photometric data for the time coverage of our Doppler maps from 
data set, see their Fig. 2c). 



Olah & Strassmeier 
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a fixed value of 300 K seems advisable. Alternatively, we 
use the minimum value increased by 500 K and an arbitrar- 
ily fixed value of 4000 K as threshold temperatures. All tem- 
perature values smaller than the respective threshold tem- 
perature provide us with a fraction value which is translated 
to spottedness in percent. The results are shown in Fig. [6] 
The threshold based on the average temperature is, as pre- 
sumed, too close to the photospheric-temperature bulge and 
sometimes exhibits large increases (around 1990 which is 
not seen in any of the other curves), induced by the scat- 
tering of the photospheric range. The minimum tempera- 
ture, increased by a specific value, seems to be valuable but 
sometimes exhibits large deviations which are not resem- 
bled in any of the other mapping lines. Overall, we conclude 
that the average temperature decreased by 300 K gives the 
most reliable threshold as it lies well between the photo- 
spheric and the spot-temperature bulge. 



Concentrating on the solid black line with filled circles 
(Fig. [6]), we notice an increase in spottedness in 1992 of 
10 percentage points in A6439 (according to a factor 2), 
also seen in the A6430 map and another small increase 
(4 to 6 percentage points) in 1995, mainly in the A6439 
map. We may compare these results with spot-filling fac- 



tors derived from TiO band modelling. Solanki & Unruh 

f04l) list filling factors of EIEri (as obtained from O'Neal 
l1. 1998 ) for the epochs 1992Mar (total 25% / minimum 
11%), 1995Jan (13%/5%) and 1995Dec (13%/5%). These 
values coincide with the increase seen in A6430 in the map 
1992 Jan as compared to 1995Dec but not to the same extent. 
However, this relation is not seen at all in A6439 and A641 1. 
Overall, no systematic changes are evident and, above all, 
no correlation with the photometric long-term trend is seen. 



3.4 Spot occurence functions 

Fig.[8]shows for all Doppler maps of the A6439 line the tem- 
perature distribution as a function of stellar longitude (left) 
and latitude (right). Starting from 11/1988 on top, each map 
is shifted by — lOOK for better viewing. Overall, we can- 
not determine any preferred surface positions of spots on 
EI Eridani (apart from the polar cap as such) which is also 
demonstrated by the "grand average Doppler map", an aver- 
age of all available Doppler maps from 1988 through 1998 
(Fig.[TT]). Does this implicate that there are no preferred lon- 
gitudes at all? This would be in contradiction to what was 
published by Berdyugina & Tuominen (1998]). Their pho- 
tometry showed active longitudes which were not synchro- 
nized with the orbital motion and were shifted by one orbital 
phase in about 2.7 years together with switching primary 
and secondary minima after 9 years (flip-flop). Fig.[9]allows 
a more detailed look at the spot occurrences as a function 
of stellar longitude. The functions are the same as in Fig. [5] 
(left) but are plotted as two mirrored contours filled in black, 
after normalization by their respective maximum temper- 
ature ("spot-occurrence function"). A large width denotes 
lower average temperature and therefore stronger spotted- 
ness. The uppermost panel shows the latitudinal locations 
of the strongest and the second strongest polar appendage 
as inspected by eye. For the years 1989 till 1994, we see a 
possible correlation of the longitudinal occurrences of po- 
lar appendages, showing a systematic drift towards higher 
longitudes, i.e. later orbital phase, indicated by solid (pri- 
mary maxima) and dotted (secondary maxima) grey lines. 
The shift amounts to 360° after about 3 years which is very 
similar to the 2.7 years found by [Berdyugina & Tuominen 



1998| ). Assuming a switch between primary and secondary 



minima to occur around 1994 (as proposed by Berdyugina 



| & Tuominen | [T998l ), we reverse the solid and dotted lines 
at 1994.4 in Fig. [9] (uppermost panel) and see that the con 
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Longitude dependence of spots (A,6439) 
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Fig. 9 The upper panel gives the longitude of the strongest (filled circles) and second strongest (if existent; open circles) 
polar appendage for each map. The tilted lines denote the suggested drifting active longitudes. Primary active longitudes are 
shown as solid lines, secondary active longitudes (on the opposite side of the star) as dotted lines. Lower panel: Longitudinal 
spot-occurrence functions. The black areas represent average temperature as a function of longitude, normalized by its 
respective maximum average temperature. A broader area denotes lower average temperature (corresponding to stronger 
spottedness) at the respective longitude. Overlapping epochs were omitted. 
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Fig. 10 Latitudinal spot-occurrence functions. The black areas mark averaged temperatures as a function of latitude. 
See Fig. ^ for further description. This display mode is useful for comparing observational results with the theoretical 
spot-probability functions presented by |Granzer et al.|pOQQ| ). 



tinuation from 1994 onwards follows roughly this pattern. 
The active longitudes noticed in Fig. [9j uppermost panel, 
are not clearly resembled by the average-temperature func- 
tions shown in the lower panels of the same figure (which 
give an average of each longitude, from —ito 90° latitude). 
Thus, we conclude that the longitudinal migration is mainly 
caused by the polar appendages. The latitudinal region of 
the polar appendages is also the most active area on EI Eri, 



as revealed by the long-term sigma maps (see Fig. 1 1 bot- 
tom). 

Fig. [To] shows the temperature averaged along latitudes. 
Clearly, the polar cap dominates the surface maps. Compare 
these results with the diagrams presented by G ranzer et al.| 
Fig. 4). They are similar only to Granzer's 0.4 



(2000 



cases, and no theoretically predicted spot-probability func- 
tions are available yet for giant stars (the results from 
Granzer et al. are for ZAIVLS stars). 



The lower panel of Fig. 1 1 shows the respective sigma 
map of a combination of all maps from 1988-1998 for the 
A6439 line. Clearly, the most and strongest variations are 
found in the appendices of the polar cap at a latitude of 60° - 
75°. Nevertheless, as demonstrated by Fig.fTOl the extent of 



tion, at ^50°. The polar cap is present in all maps|^If the 
polar cap would vanish for a certain epoch, we would ex- 
pect the respective line profiles to be more pointed, i.e. less 
shallow. While all averaged profiles are very homogeneous, 
a few deviate slightly in line depth. This behaviour is, how- 
ever, not followed in all three mapping lines simultaneously 
and is therefore likely due to noise. In fact, none of the ex- 
tracted parameters exhibits a clear systematic variation in 
all three mapping lines. 

4 Discussion 

Our long-term DI study revealed that the brightness vari- 
ations (exhibiting a period of roughly 12 years) are not 
resembled by any of the parameters extracted from our 
Doppler maps nor by the directly measured line equivalent 
widths. This is in contrast to the solar analog where the to- 
tal spottedness varies with the activity cycle. Thus, we are 
asked to explain what exactly causes the long-term varia- 
tions of the mean brightness and why are these variations 
not reflected in our Doppler maps. 



the polar cap remains very stable and ends, without excep- parameters 



9 I 

Unlike wrongly claimed from preliminary results by Washuettl & 

Strassmeier^ pOQ2j due to an ill-posed Doppler map and inconsistent DI 



www.an-journal.org 



© WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



8 



A. Washuettl et al.: The active binary star EI Eridani: Long-term Doppler imaging 



50 

40 

30 

20 

10 


50 

40 

' 30 

20 

10 


50 

40 

30 

20 

10 



Fractional spottedness for various threshold values 
1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 



average temp. - 300K [%\ 
Linear fit to av. temp. - 300K 
o - ■ -o average temp. [%] 
• minimum temp. + 500K [%] 
4^-4 4000K[%] 

average temp./lOO [K] 





A 4 



48000 49000 50000 51000 
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Fig. 6 Fractional spottedness in % of the entire sphere for 
all Doppler maps as a function of time and various thresh- 
hold values (explained in the insert). The solid black line 
with filled circles, i.e. the average temperature minus 300 K 
(solid black line), gives the most consistent results. The grey 
line denotes the correponding linear fit, see Tab. [6] For com- 
parison, the average temperatur (divided by 100) is plotted. 



4.1 Short-term variations 

The photometric amplitude varies such that at the begin of 
the photometric minima in 1984/85, as well as in 1995/96, 
the amplitude was smallest and turned largest shortly there- 
after. A clear correlation with the long-term trend is not ob- 

left up- 



vious (see Fig. 2c in Olah & Strassmeier 2002 



per panel). These short-term brightness variations are as- 
sumed to be caused by the same spots moving across the 
visible disk of the star. However, differences in the ampli- 
tude of those variations do not necessarily reflect a larger or 
a smaller number of spots or spottedness. They can also be 
a sign of a homogeneous or otherwise a concentrated spot 
distribution on the star's surface. However, the short-term 
brightness changes, if caused by spots, should be directly 
related to the longitude functions in Fig.[9j and should there- 
fore be seen. Most likely, the brightness changes are only 
caused by those parts of the stellar surface that are not cir- 



Temperature as a function of longitude Temperature as a function of latitude 




Fig. 8 Temporal variations of surface temperature as a 
function of longitude (left panel) and latitude (right panel). 
Each function is normalized by the average of all functions 
and shifted by 0.1 for better viewing. To the right, year and 
month of observation is unscrambled. The overall temper- 
ature profile shows only small and uniform variations. The 
largest variations are found at high latitudes (65-80°), as 



also seen in Fig. 1 1 (bottom panel). 



Average map 1989-1998 (A6439) 




El Er 
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Sigma map 1989-1998 (A6439) 
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Fig. 11 Top: the grand average Doppler map for 1988- 
1998 (A6439). No preferred longitudes in the orbital refer- 
ence frame are seen. Bottom: Residual map for 1988-1998 
(A6439). 



© WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN () 



9 



cumpolar. Therefore, we recreated the plot displaying the 
longitudinal spot-occurrence functions with the high lati- 
tudes (90° — i) excluded. However, a relationship with the 
brightness variation is still not obvious. After all, for epochs 
with a small number of photometric data points, the disad- 
vantageous rotational period of almost two days may give 
an explanation: most photometric data of EI Eri were col- 
lected by automatic telescopes which normally execute the 
same targets every night. This brings about a similar time 
of observation of a specific target for each night and, in the 
case of EI Eridani, a similar orbital epoch which will, for a 
small data sample, mimic smaller photometric amplitudes. 

4.2 Long-term variations 

A decrease or increase in total brightness - if caused by 
a larger spottedness or, like on the Sun, by a higher facu- 
lae emission - should accordingly be reflected by the spot- 
tedness as extracted from the Doppler maps, and thereby 
also by the average temperature and the equivalent width. 
However, none of it is the case with EI Eri. One possi- 
ble explanation is that small inaccuracies in the normaliza- 
tion of the spectra (continuum fitting) could potentially in- 
fluence the temperature span and thereby introduce noise. 
Most likely though is that the cyclic brightness variation 
is caused by a fluctuation of small unresolved spots while 
the fractional coverage of large-scale spots remains basi- 
cally constant. Doppler images do a good job in catching 
starspots that significantly modulate the line profile, but 
it is extremely difficult to detect a background of small 
starspots evenly distributed over the stellar surface. |Jef-| 
fers et a l.] (^2006) obtained spectrophotometric eclipse light 
curves with HST/STIS and mapped the primary F9 compo- 
nent of the F9V+K4V-IV binary S V Cam. They found that 
the observed surface flux from the eclipsed low-latitude re- 
gions of the F9 primary was 30% lower than what is pre- 
dicted from model atmospheres and the Hipparcos paral- 
lax. This could be explained only if there is about a third of 
the eclipsed region covered with unresolved cool starspots. 
Even if this "dark" spot component is assumed to also ex- 
ist over the rest of the surface, a huge cap-like polar spot 
down to a latitude of 48° was required in order to fit the 
eclipse light curves. [Solanki & Unruh| ( |2004[ ) found that 
even for heavily spotted (hypothetical) stars, a large frac- 
tion of the spots are smaller than the current resolution limit 
of Doppler images. Thus, it is possible that a substantial 
part of the spot coverage is not resembled by our measure- 
ments of fractional spottedness (as presented in Fig.|6]). Spot 
covering fractions can also be derived from TiO band mod- 
elling, a method that does not suffer from the resolution 
limitation. It would therefore be essential to monitor and 
compare fractional spottedness as derived from DI with TiO 
filling factors. However, |Berdyugina| (|2002f) obtained filling 
factors for the RS CVn system II Peg that agree with the 
corresponding Doppler images without requiring additional 
(unresolved) starspots. As a final alternative, we could only 
assume that the long-term variation of the total mean bright- 



ness is not caused by the degree of spottedness. In this case, 
it would be questionable in what way the photometric cycle 
corresponds to a magnetic cycle. 

4.3 Comparison with the Sun 

On the Sun, the solar cycle is visible not only in the vi- 



sual brightness (see e.g. |Willson & Mordvinov||2003j but 
also in many other tracers including total irradiance, radio 
and total magnetic flux. Call H & K emission, in mag- 
netic features like sunspots, and others (see Harvey & White] 



1999). The sunspot number, though, is the most prominent 



one. Most indices vary roughly in phase with the sunspot 
number. However, the modulation amplitude varies widely 
between different indices. In visible light, the amplitude is 
marginal and vague, and most likely not obvious to a casual 
observer. 

The sunspot maximum correlates with a maximum in 
total irradiance/brightness, Ca II H & K core emission and 
total magnetic flux ( Harvey & White||1999[ [Radick et aL] 



1998). Thus, the Sun is brighter during sunspot maximum 



which is caused by the dominating emission from faculae. 
This behaviour is also found on other old stars ( [Radick et al.| 
|1998| ). Both, the Sun and EI Eri, tend to become bluer as 
they get brighter, which suggests that EI Eri might exhibit a 
higher level of faculae emission during activity (i.e. photo- 
metric) maximum. On the Sun, however, the emergence of 
faculae is usually associated with spots. 

A model presented by Solanki & Unruh ( [1998 ) indicates 
that the radiative properties of magnetic features on the so- 
lar surface (i.e. faculae and sunspots) provide the dominant 
contribution to irradiance variations on a solar-cycle time- 
scale, suggesting that such features are also responsible for 
the visual brightness variations on EI Eri. However, the to- 
tal irradiance variations seem to be strongly dependent on 
the inclination angle of the star and, for the Sun, could in- 
crease by a factor exceeding 6 when viewed from a heli- 
ographic latitude of 60° (see Radick et al.|[T998 ). Possibly, 
the variations in total spottedness on EI Eri are minute and 
not resolved by our Doppler images while the accompany- 
ing faculae could have an enhanced effect on the brightness 
and photometric colour. 

4.4 Comparison with other stars 

The surprise that the proposed photometric cycle does not 
have a corresponding spot or magnetic cycle might get sup- 
port by results from other stars. Donati et al.| ( [2003 j ) pre- 
sented long-term magnetic surface images for the young K 
dwarfs AB Doradus and LQ Hydrae and for the RS CVn- 
type Kl subgiant HR 1099 and found that large regions with 
predominantly azimuthal magnetic fields are continuously 
present at the surface of these stars. Long-term structural 
changes take place but reflect no more than the limited life- 
time of the corresponding surface structures. An active lon- 
gitude is not confirmed on HR 1099. Furthermore, no clear 
secular change is detected in the axisymmetric component 
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of the magnetic field and in particular, no global polarity 
switch was observed in the field ring pattern of HR 1099, 



nor in the high-latitude azimuthal-field ring of AB Dor. Petit 



etal.| ([2Q04 ) report that, on HR 1099, the small-scale bright- 
ness and magnetic patterns undergo major changes within a 
time-scale of 4-6 weeks, while the largest structures remain 
stable over several years. 

Another case, the K2 IV RS CVn star II Peg, exhibits - 
like EI Eri - migrating active longitudes including the flip- 
flop phenomenon. However, the total spot area is approxi- 
mately constant during this cycle. This means that II Peg's 
activity cycle is only expressed as the spot area evolution 
within the active longitudes, i.e. as a rearrangement of the 
otherwise nearly constant amount of the spot area ( ,Berdyug-, 
ina etal.|1999l ). 



Vogt et al. (19991, having observed HR1099 for 11 



years, suggest starspots to be stellar analogs of solar coronal 
hole structures and expect a dynamo cycle to be manifested 
by periodicity in the area of the polar spot - an expectation 
not confirmed by EI Eridani. 

4.5 Preferred longitudes 

As pointed out above, EI Eri seems to exhibit preferred lon- 
gitudes that migrate within the orbital reference frame, sup- 
porting the finding from |Berdyugina & Tuominen| ( |1998| ). 
The time it takes for an active longitude to shift by one 
orbital phase (360°) is about 3 years which is close to the 
short-period cycle found bylQlah & Strassmeier ( 2002 |). Can 
this shift be caused by differential rotation? A evaluation 
of differential rotation on EI Eri suggests a value of a ^ 



0.037 (i.e. solar, see the forthcoming paper Washuettl et al. 
2009| ). If this value is real, the equatorial spots overlap the 
polar region in about 50 days. Particularly, differential rota- 
tion then cannot account for the longitudinal migration of 
spots as suggested by the phase drift of the migrating pho- 
tometric wave minimum - a result also found on HR 1099 
( |Vogt et al.|1999 ). During the time the equator overlaps the 
pole (^50 days), the preferred latitudes shift only by about 
16°. After one orbital rotation, this shift amounts to ^ 0.65°, 
while the shift of the equatorial region comes to about 14° 
during the same time. Therefore, considering both differ- 
ential rotation and the existence of stable, preferred longi- 
tudes, the spots constituting the high-latitude appendages 
must either be very short lived (otherwise they would be 
moved out of the preferred longitudes by the slowerer ro- 
tating polar region) or they are anchored in a deeper layer 
that is not bound to the differentially rotating surface. Possi- 
bly, the appendages at the preferred longitudes (or perhaps 
all spots in the polar region) are the footprints of a strong 
dipole magnetic field that is firmly anchored to the synchro- 
nized deeper core of the star. The magnetic energy density 
of such a dipole could be stronger than the kinetic energy 
density of the differential shearing motions of the gas. This 
emergent dipole field would then dominate the gas motions 
at high latitudes, producing apparent synchronization of the 
high-latitude features with the orbit. 



Active longitudes seem to be a common feature on ac- 
tive close-binary stars. Such a phenomenon has been no- 
ticed using the technique of Doppler imaging on, e.g., the 
active giant DM UMa (Hatzes| [T995|), on the two main- 
sequence components of CrB (Strassmeier & Rice 



2003 ), on the two pre-main sequence components of V824 
Ara jHatzes & Kurster| |1999i Strassmeier & Rice 7000] ) 



and the active giant UZ Lib ( Q lah & Strassmeier 2002| ), as 
well as by long-term photometric analysis on, e.g., RT Lac 
( [gakirli et al.pQQ3] ), UZ Lib (Qlah et al.|2Q0 2a) and on the 



RS CVn stars II Peg, a Gem and HR 7275jBerdyugina& 
|Tuominen|1998| ). [Holzwarth & Schussler|p003| presented 



a numerical investigation for fast-rotating solar-like binaries 
which showed that, although the magnitude of tidal effects 
is small, they nevertheless lead to the formation of clusters 
of flux tube eruptions at preferred longitudes on opposite 
hemispheres and synchronized with the orbital motion. This 



finding is supported by simplified calculations by Moss & 



Tuominen ( 1997 ) who showed that synchronized close late- 
type binaries can be expected to exhibit large-scale non- 
axisymmetric magnetic fields with maxima at the longitudes 
corresponding to the two conjunctions (again locked within 
the orbital reference frame). None of these investigations 
can reproduce migrating active longitudes and the flip-flop 
phenomenon as seen on EI Eri. For a possible theoretical 
explanation of the flip-flop phenomenon on the single giant 
FK Com see |Elstner & Korhonen| ( |2005] ) (see also |Korhonen| 
|& Elstner|2005l [Tuominen et al.|2002] ). 



5 Summary 

The key results from the Doppler-imaging analysis are as 
follows: 

- The general morphology of the spot pattern persisted for 
more than 10 years. 

- The polar spot (AT ^ 1500 K) is stable. No decay and 
reemergence due to, e.g., a polarity inversion was seen. 
However, it changes its shape on short time-scales (one 
week). The size of the polar cap is stable and reaches, 
without exception, down to a latitude of ^50°. 

- The most active surface region is the latitude between 
60° and 75°; the polar-spot appendages in this region 
seem to be responsible for the light variability and espe- 
cially for the migrating active longitudes. 

- Low latitude spots occur and decay on short time- scales 
(less than a week) and are much less pronounced (AT ^ 
500 K). 

- The polar appendages seem to appear on preferred lon- 
gitudes and drift towards larger longitudes with respect 
to the apsidal line of the binary system with a period 
of ^3 yr, thereby confirming the migrating active longi- 
tudes found earlier by |Berdyugina & Tuominen] ([1998 ). 
No fixed preferred longitudes locked to the orbital ref- 
erence frame of the system are seen though. 

- No correlation is seen between the photometric long- 
term cycle and any parameter extracted from our spectra 



© WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN () 



11 



and Doppler maps; this includes equivalent width, lat- 
itude and longitude function, average temperature and 
spottedness. 
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